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Abstract: We investigate the possibility that the process of (j) photoproduction may have 
a significant hard perturbative QCD component. This suggestion is based on a study of the 
energy dependence of the forward photoproduction cross section followed by a calculation 
where we show that a coherent sum of the pQCD and conventional soft Pomeron contribu- 
tions provides an excellent reproduction of the experimental data. Our results suggest that 
the transition from the predominantly soft photoproduction of light p and u vector mesons 
to the predominantly hard photoproduction of heavy J/\I/ and T is smooth and gradual, 
similar to the transition observed in deep inelastic scattering studies of the proton structure 
function in the small x limit. Our predictions for higher HERA energies are presented. 
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1 Introduction 



Over the past few years we have seen a vigorous phenomenological investigation of the 
Pomeron through the study of hadronic total, elastic and diffractive cross sections, as well 
as the study of the proton deep inelastic scattering (DIS) structure function. In particu- 
lar, Donnachie and Landshoff (DL) have promoted an appealing and very simple Regge 
parameterization of the total hadronic cross sections in which 



The two key ingredients of this approach are the Regge trajectories 

a R (t) = a R (0) + a' R t, (2) 
where or(0) = 1 — r\ and the Pomeron trajectory, which dominates at high energies, 

a P (t) = a P (0) + a' P t } (3) 

where ocp(0) = 1 + A. The DL study establishes universal values A = 0.0808 and rj = 0.4525. 
This study is supplemented by the analysis of Block, Kang and White |2| who determine the 
slope a' P of the Pomeron trajectory to be a' P = 0.2 GeV -2 . 

The same approach is also applicable to the analysis of real photoproduction and of the 
proton DIS structure function 0. While the energy dependence of the photoproduction total 
cross section follows the DL pattern, it has been observed [§]] that F 2 (x,Q 2 ) behaves, for 
small enough x, like x~ x , where A is slowly growing with Q 2 . The growth of A is associated 
with the behavior of the gluon structure function in the small x limit 

din (xG(x,Q 2 )) 

dln(l/x) ' 1 ' 



It has been recognized for quite some time that the transition from the predominantly 
soft real photoproduction (Q 2 = 0) to the predominantly hard DIS processes, with high 
enough Q 2 , is smooth and gradual [§, §, ||, [/J. This observation, regardless of its theoretical 
interpretation, is evident once we examine the energy dependence of F 2 (x, Q 2 ) in the small x 
limit with Q 2 ranging from zero to a few GeV 2 . It is also well known that real photoproduc- 
tion of light vector mesons, p and u, is dominated by a soft Pomeron exchange [§ 0, whereas 
photoproduction of heavy vectors, J and T, is well reproduced by a perturbative QCD 
(pQCD) calculation [□], [12], |13|, where My/ '4 replaces Q 2 as a measure of the process 
hardness. Although M 2 /4 is a discrete variable, while Q 2 is a continuous DIS variable, it 
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is interesting to check if the transition from soft to hard photoproduction of vector mesons 
follows the behavior pattern observed in DIS. To this end, the study of photoproduction is 
particularly instructive, as Mj/4 = 0.26 GeV 2 , while we know that the energy dependence 
of F 2 , with small x and Q 2 as low as 0.2 — 0.3 GeV 2 , is steeper than the energy dependence 
of cr tot (7p). 

Our investigation is susceptible to both experimental and theoretical uncertainties. Ex- 
perimentally, a systematic study of the integrated photoproduction cross section and the 
forward differential cross section slope [[14], [15], [16], |T^| are not very reliable as different exper- 
imental groups have utilized different, not always mutually consistent, methods to extract 
and relate these quantities. To overcome this difficulty, we have analyzed the measured 
differential cross sections rather than integrated quantities. Even so, the two higher energy 
data points fT6| , |TT| are averaged over wide energy bins. This, combined with the over- 
all poor quality of the reported data, may make a detailed analysis non conclusive at this 
stage. Theoretically, since we wish to utilize the gluon structure function at low Q 2 , the 
calculation of the hard component requires some clarifications. Technically, a pQCD calcu- 
lation of [da^p — > 4>p)/dt] t= o depends on our knowledge of the gluon structure function at 
Q 2 = 0.26 GeV 2 . Such information requires an extrapolation of a given parton distribution 
below its initial evolution threshold Qq. For this purpose we adopt a linear extrapolation 
which was successfully utilized in previous calculations P, [H|. As we shall see, there is a 



significant difference between the MRST jl9|] and GRV98 f2(| input gluon distributions. We 



have chosen to use the GRV98 distribution and shall explain our motivations for doing so. 

The purpose of this letter is to examine these issues in some detail from different points 
of view. We present an analysis of the existing photoproduction forward differential cross 
section data which suggests an energy dependence which is steeper than the typical energy 
dependence associated with the soft Pomeron [l|]. We then present a pQCD calculation 
from which we deduce that the hard component is responsible for about a quarter of the 
photoproduction amplitude in the forward direction at presently available energies. We then 
proceed to show that a coherent sum of the calculated pQCD amplitude and a conventional 
soft Pomeron exchange contribution provides an excellent reproduction of the available data 



14, 15, 16, 17 



2 Data analysis 

Our data analysis investigates whether the photoproduction cross section follows a power 
dependence on the cm. energy W, and whether this power is larger than the value deter- 
mined from the energy dependence of the total cross section. Following DL [|]] and Block et 
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al. ||, we expect the (ft photoproduction cross section to behave like PU 4A and the forward 
slope to behave like 4a' P In W. The analysis of (ft photoproduction data is seemingly easy, 
as this process proceeds exclusively through Pomeron exchange, since the various Regge 
exchanges cancel each other. The problem is that the published analysis M |T7[ depends 
on a comparison between integrated cross section data taken by different groups who have 
utilized different, and not always mutually consistent, procedures. In addition, because the 
(ft forward slope is shrinking, the interpretation of the integrated cross section behaving as a 
fixed power of W is somewhat ambiguous. In order to bypass these difficulties, we have lim- 
ited ourselves to the analysis of the individual differential cross sections da/dt as reported by 



the experimental groups |L4], [L^, [16], [T7|]. We have used data with W > 6 GeV corresponding 
to x < 0.025. 

Fig. 1 shows that (da/dt)o in the available energy range is, indeed, well fitted by an 
effective power of W, 



/ da(-fp -¥ (ftp) 



AW AX . (5) 



=o 



Our best fit, for 5 data points with 6.7 < W < 70 GeV, has an excellent x 2 jn.d.j. = 0.22, 
corresponding to A = 0.76 ± 0.09 and A = 0.135 ± 0.012. For comparison we show also a 
fit where we fix the power to its DL value 4A = 0.3232 and obtain A = 1.21 ± 0.06 with 
X 2 /n.d.f. = 0.92. This is a lesser quality fit, but it cannot be discredited. Our inability to 
determine the power unambiguously results from the big error coupled to the ZEUS high 



energy data point [17]. Even so, it is clear that an improvement of the HERA data point 
at < W >= 70 GeV and relevant data at higher energies will enable us to conclusively 
distinguish between a DL type interpretation and ours. 

In order to further examine the suggestion that the dependence of (ft photoproduction on 
W is steeper than the behavior implied by a soft Pomeron exchange, we have studied two 
(related) ratios 

_ (da/dt) 
(dcr/dt) 

where (Jpipfp) is the soft Pomeron exchange contribution to cr tot (7p). Both ratios are shown 
in Fig. 2. A careful study of these ratios is of interest as both crp^p) and a 2 ot ((ftp) behave 
as W 4A . Using the DL parameterization |IJ we have 

0.1616 



and 



a P (jp) = 67.7 j fib 
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and, with aid of the additive quark model relation between cj)p, Kp and irp cross sections, 



, w . 0.1616 , w , -0.4525 

*«(*p) = 10.01 (-) +1.51 (-) mb, (9) 

with Wq = 1 GeV. A predominantly soft production mechanism for 'jp — > <ftp means that 
both ratios presented in Fig. 2 are constants. The best fits that we have obtained imply 
that Ri behaves as w°- 27S±om6 and R 2 behaves as w°- 215±0 - 020 . Constant ratios provide 
marginally acceptable fits to R\ and R 2 -, but this option cannot be definitely excluded. 

To summarize: The available experimental data on [da{^fp — > <fip) / dt) Q is consistently well 
described as having an energy dependence steeper than the one implied by a soft Pomeron 
exchange. Nevertheless, the assumption of a pure soft production mechanism cannot be 
unambiguously eliminated. 



3 A pQCD calculation 



Our pQCD calculation of the forward photoproduction follows earlier pQCD calculations 
of the forward photoproduction cross section of heavy vector mesons [IIJ [11], [12], |I3|| . These 
calculations are considerably simplified once we assume a non-relativistic wave function for 
those vector meson states. This assumption, which is also valid for 0, enables us to write a 
leading-order expression 



da(jp —> (ftp) 
dt 



t=o 



xG x, 



Ml 



(10) 



where x = (M^/W) 2 and Tf e is the partial decay width of (p — ► e + e~ . In the following we 
follow Refs. Jl2|][T3|] and, after calculating the cross section resulting from the imaginary 
forward amplitude in leading order, we correct for the real part of the amplitude and for 
higher orders. 

The basic problem with our suggested calculation is that we depend on the knowledge 
of the gluon structure function xG(x, Q 2 ) at Q 2 = M|/4 = 0.26 GeV 2 . This Q 2 value 
is well below Q^, the initial evolution threshold, used in the updated parton distribution 
parameterizations [|T^][^D| which took into account the behavior of F 2 and its logarithmic 
derivatives at small Q 2 . We recall the general property of the gluon structure function which 
is linear in Q 2 in the limit of very small virtuality and use, accordingly, a linear extrapolation 



Q2 

xG(x,Q 2 ) = -—^xG(x,Ql) 



Q 2 <Ql 
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Figure 2: R\ and R2 data and best fits. 
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This approximation has been successfully used in previous theoretical DIS investigations 
[0 10 m which the knowledge of xG in the small Q 2 region was required. 



The results of our gluon structure function extrapolation for MRST |T^[ and GRV98 |f2D 
at the relevant Q 2 = 0.26 GeV 2 value, are shown in Fig. 3. Note that Q 2 = 1.2 GeV 2 for 
MRST and 0.8 GeV 2 for GRV98. Fig. 3 shows also the Pomeron term of ALLM97 ||. In 
the following we have used the GRV98 extrapolated distribution. Our motivation is double 
folded: 

1) From a practical point of view, the MRST gluon structure function is considerably smaller 
than GRV98. If we adopt the MRST distribution, we get a diminishing small hard contribu- 
tion for photoproduction in the W range of interest. If our data analysis is substantiated, 
we would then be left with the need for some explanation for the observed energy dependence 
of the data. 

2) Theoretically, we note that MRST is very close to ALLM97 in the x region of interest. 
Our interpretation is that the MRST input at small Q 2 is predominantly soft, such as is 
ALLM97, and thus, even though perfectly legitimate, it is less suitable for our analysis. 



We follow Ref . [|T2"|, |13|| and consider the following corrections to the leading order cross 
section written in Eq. (|TT)|): 

1) Eq. ( |10"D corresponds to the imaginary part of the forward amplitude and should be cor- 
rected by (1 + p 2 ), where p = KeA/lmA, and A is the amplitude under consideration. We 
recall that p ~ 7rA/2, where, in our case, A is given by Eq. ([|). We note that A = 0.145 for 
W = 70 GeV and decreases very slowly as the energy decreases. 

2) The next to leading order corrections are estimated by [1 + 0.5as{M^/4:)]. 

3) Relativistic corrections and the effects of intermediate off diagonal partons in <p photo- 
production are rather small and have been neglected. 

With these effects put together, our overall correction factor is C = 1.21 at the high 
energy end of our data, decreasing slowly to 1.16 in the lowest energy. The energy dependence 
of our calculated hard cross section is given in Table [I] and displayed in Fig. 4. 



4 A hybrid pQCD and soft Pomeron model 

Although the pQCD contributions calculated in the previous section provide a significant 
contribution to photoproduction, there is no doubt that, in the energy range under con- 
sideration, the leading production mechanism is a soft Pomeron exchange. Accordingly, we 
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Figure 3: ALLM97 and the extrapolated MRST and GRV98 parameterizations for the gluon 
structure function xG(x, Q 2 ) at Q 2 = 0.26 GeV 2 . 
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\dtJo 
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3.46 


0.61 
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3.80 


0.67 


250 




8.20 


4.08 


0.71 



Table 1: {dcr/dt) t=0 data and calculations. 




attempt to fit the data with a simple hybrid two- component model with the following pre- 
scriptions: 

1) The first component is a soft DL Pomeron with an ap(t) intercept of 0.0808, namely 

Soft w 0.3232 

' <12) 

/ j \ Hard 

2) A hard pQCD component > as calculated in the previous section. 

3) A Coherent sum of the two component amplitudes. 

We fit the 5 data points of (da/dt)o with one parameter, the normalization As of the DL 
Pomeron. We obtain a best fit value A s = [0.83 ± 0.02] (/i6) 5 with x 2 /n.d.f. = 0.54. Our 
fitted cross sections (called FM) are presented in Table [I] and displayed in Fig. 4, both table 
and figure showing also the pQCD and soft cross sections and our predictions for higher 
HERA energies. Our fit should be compared with the power fit, presented earlier, which 
has a x 2 / n -d-f- = 0.22 and the conventional DL fit which has \ 2 /n.d.f. = 0.92. Clearly, 
the presently available data is not sufficient to rule out any of these options. Once again 
we note that this ambiguity results from the big error associated to the < W >= 70 GeV 
point. Improvement of the quality of this point and additional HERA data will enable a 
more discriminative analysis. 

The hybrid model that we have just suggested can be further examined by considering 



the differential cross sections. Such data is available 16, 17 at < W >= 13.7 and 70 GeV. 
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For the purpose of our calculations we need to know the t-dependence of both the pQCD 
and soft Pomeron amplitudes. To this end we assume each of these dependences to be 
approximated by an exponential, 



Fi = Fi(0)e^* i = S,H. (13) 

This simple approximation is sufficient here, considering the quality of the available data on 
photoproduction. We then proceed as follows: 

1) We take for the pQCD amplitude an energy independent exponential slope. We derive 
its value from the high energy differential cross section on J/ty photoproduction combined 
with the observation pi] that this slope is energy independent, corresponding to a flat 
hard Pomeron. In our analysis we have taken Bh = 4.6GeV~ 2 , which corresponds to the 
HI measurement p2f . An equally good fit is obtained also with the ZEUS value [^3| of 
Bh = 4.0 GeV -2 . We have also treated Bh as a free parameter. 

2) For the soft Pomeron amplitude we assume a conventional Regge type exponential slope 
depending on two fitted parameters, 

/ W\ 2 

B s =B* + 2a' P \n — . (14) 



W 

The Bq approximate value is known from the phenomenology of (ft photoproduction at low 
energies and a' P is approximately known from high energy hadron-hadron phenomenology 
M. As we show below, the fitted values of these parameters are in excellent agreement with 
our expectations. 

This hybrid model was fitted to reproduce 13 (da/dt) data points measured at < W >= 
13.7 and 70 GeV. We obtain B$ = 4.20 ±0.90 GeV" 2 and a' P = 0.20 ±0.08 with ^jn.d.f. = 
0.65. Our results are in an excellent agreement with soft hadron Regge analysis 0. Our soft 
exponential slopes extrapolate well into the predominantly soft low energy experimental data 
on <f) phtoproduction The data and our fit are shown in Fig. 5. A slight improvement 
of the fit is obtained if we take Bh as a free parameter and obtain Bh = 3.50 ± 0.60 GeV -2 , 
with the other parameters changing insignificantly. 

Although our (da/dt) fit corroborates our proposition that high energy <ft real photopro- 
duction has a significant hard component, we do not consider our success to be decisive. The 
< W >= 70 GeV data, on its own, can be equally well described by a single DL type soft 
component. Our two component model is appropriate to describe the < W >= 13.7 GeV 
data including t > 0.4 GeV 2 , where a single exponent is not sufficient. However, we caution 
against reaching too strong a conclusion from a single measurement. Clearly, and not only 
for the purpose of our analysis, additional knowledge on higher t behavior will help to clarify 
the picture. 
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Figure 5: Data and our calculations for (da/dt) at < W >= 13.7 and 70 GeV. 
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5 Discussion 



Following are the main conclusions of our study and some general remarks: 

1) The data analysis of forward real <fr photoproduction suggests the existence of a significant 
hard pQCD contribution in the energy range of 6 < W < 70 GeV. 

2) The above suggestion is corroborated by a pQCD calculation of [da{^p — > <j)p)/dt]o using 
the GRV98 gluon distribution extrapolated to Q 2 = 0.26 GeV 2 . 

3) A hybrid model, in which we coherently add the parameter free pQCD hard component 
and a DL type soft component, provides an excellent overall reproduction of the data. We 
note that the ratio of hard to soft contribution in our model is compatible with the ratio 
obtained in a detailed study of DIS in comparable small Q 2 values. 

4) Our hybrid model is significantly different from the two Pomeron model suggested in Ref. 
[0. In our model the t = intercept of the hard effective Pomeron is given by Eq. (|]) and 
as such it is a moving pole. In the model of Ref. |7], the hard Pomeron is a fixed pole with 
a comparatively high A. 

5) Theoretically, the validity of our calculation rests on (i) the legitimacy of our Q 2 extrap- 
olation of xG(x, Q 2 ) below Ql ; and (ii) our choice of GRV98 for the input gluon structure 
function. 

6) Our overall analysis strongly supports the existence of a significant hard component con- 
tributing to 7p — > (pp. However, a decisive quantitative conclusion depends on improving 
and extending the HERA data on <fi photoproduction. 
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